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Abstract 

We study the Yang-Lee zeros of a random matrix partition function with the global 
symmetries of the QCD partition function. We consider both zeros in the complex chem- 
ical potential plane and in the complex mass plane. In both cases we find that the zeros 
are located on a curve. In the thermodynamic limit, the zeros appear to merge to form 
a cut. The shape of this limiting curve can be obtained from a saddle-point analysis of 
the partition function. An explicit solution for the line of zeros in the complex chemical 
potential plane at zero mass is given in the form of a transcendental equation. 
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1. Since the work of Yang and Lee |l|], zeros of the partition function have become an 
important tool for the study of phase transitions. (See, e.g., p], [3], |J].) Most notably, Yang 
and Lee proved the theorem that the thermodynamic limit of the free energy is analytic 
in any region of the complex fugacity plane which contains no zeros. Below, we use the 
term Yang-Lee zeros as a generic name for zeros of the partition function. (Note, however, 
that zeros in the complex temperature plane were first discussed in 0, |(J.) Historically, 
Yang-Lee zeros have been studied primarily within the framework of statistical mechanics. 
Yang-Lee zeros of the Ising model have received particular attention. (See p] for recent 
work on this topic.) Recently, however, Yang-Lee zeros were studied in the context of the 



lattice QCD (QED) partition function for non-zero chemical potential ||, |9], [ID], [11], [12 
Because of the difficulty of simulating the partition function of lattice QCD for non- 
zero chemical potential fl3| , these studies are still in an exploratory stage. In this work 



we study Yang-Lee zeros for a much simpler random matrix partition function which 
possesses the global symmetries of the QCD partition function. This enables us to obtain 
the Yang-Lee zeros both numerically and analytically. 

The continuum Euclidean QCD partition function for a quark mass matrix, m, and 
chemical potential fi can be written as 

Z(m, n) = (det(7 • D + m + i^o))s QGD » CO 

where 7 • D is the Euclidean Dirac operator and 7 M are the Euclidean Dirac matrices. 
The average is over the Yang-Mills action. In lattice QCD, the chemical potential is 
incorporated by including a factor e M for links forward in time and a factor e~ M for links 
backward in time. Since the lattice QCD Dirac operator is a finite matrix, the partition 
function is a polynomial in m and e M . This makes it possible to study zeros of the partition 
function in the complex fugacity plane. The difficulty with lattice studies is that the 
coefficients of this polynomial are necessarily obtained numerically. This is problematic 
because the zeros of a high-order polynomial are notoriously sensitive to the values of its 
coefficients. 

As is well known, the zeros of Z in m are closely related to the chiral order parameter 
||. If these zeros are located at z^, the partition function ([]]) can be written as 

Z{m) = H{m - z k ) . (2) 

k 



Therefore, the chiral condensate (for one flavor) is given as 

EW.i^logZW^i:^ (3) 

We consider E(m) as a function of complex m. Provided that these zeros merge to form 
a cut in the thermodynamic limit, S(m) will show a discontinuity each time m crosses 
such a line. Similarly, the 'baryon density' defined as ns = <9 M logZ shows a discontinuity 
if [i crosses a line of zeros in the complex /i-plane. 

2. In this work we study the Yang-Lee zeros for a random matrix model of the QCD 
partition function. In essence, we replace the matrix elements of the Dirac operator by 
Gaussian distributed random variables consistent with the global symmetries of the QCD 
partition function (see jnj for a review). For Nf flavors this partition function is defined 



as 

Nf 



Z {m ,,)^jDCP{Cm^[j' +ii iC m +") , (4) 

where C is an arbitrary complex n x n matrix and DC the Haar measure. The probability 
distribution P{C) is given by 

P{C) = exp(-nS 2 TrCC t ) . (5) 

We emphasize that this partition function is a schematic model of the lattice QCD par- 
tition function. For example, it does not reproduce the zero temperature Fermi-Dirac 
distribution. However, this partition function shares certain important features with the 
QCD partition function. We mention four properties: 1. This model shows a chiral phase 
transition as a function of the chemical potential. Chiral symmetry is broken at zero 
chemical potential and is restored above some critical fi. At fi = the chiral condensate 
is E. 2. The Dirac operator is non-hermitean with eigenvalues distributed in the complex 



plane. 3. As first shown by Stephanov |15] and confirmed in [16], the quenched limit of 
this partition function is obtained as the limit Nf — > of (|j) with the fermion determinant 
replaced by its absolute value. 4. If /i = inT, it can be interpreted as a model of the 
QCD partition function with only the lowest Matsubara frequency included |TT| , |T8| , [T9| . 
This model shows a second order phase transition at 7rTS = 1. 

Obviously, the partition function (|4|) is a polynomial in m and fi. The advantage of 
studying Yang-Lee zeros in this model is that the coefficients of the polynomial can be 



obtained analytically. This eliminates a major source of error in the calculation. Here, 
we restrict our attention to a single flavor, and we adopt units in which the parameter E 
is equal to 1. 



Using standard methods [17|, [2(|, this partition function can be rewritten as 



Z N {m^) ~ J dada*e- n ^det n ( ° + ™ ■ (6) 

In the thermodynamic limit this partition function can be evaluated by a saddle point 
approximation. We obtain a non-trivial saddle point 

a = a, 

er(m + a) 2 — fi 2 a = m + a . (7) 

For > 0.527- ■ -, the trivial saddle point, a = 0, is dominant for m = [IS]. Chiral 
symmetry is restored above this point. The point at which the phase change occurs is 
determined by the solution to a transcendental equation. For complex fi it is given by 

Re[/i 2 + log(/i 2 )] = -l . (8) 

We wish to note two special solutions of this equation, /i = 0.527 • • • discussed above and 
jM — i. The latter solution corresponds to the critical temperature nT = 1 found in |T7 



We will show below that the zeros in fi of Z(m = 0,/x) for finite n lie along this curve. 
The discriminant of the cubic equation (0) is given by 

D 3 = l(mV " ™W - 5/, 2 + + /x 2 ) 3 ) • (9) 

We will also show that the endpoints of a line of Yang-Lee zeros are given by the points 
where this discriminant vanishes. 

The partition function (Bp can be written as a polynomial in fi and m by a straight- 
forward expansion of the determinant in powers of \x and m. With the help of a number 
of combinatorial identities, we arrive at 

„. . ttJV! "" WHVy (N-j)\ 

z(m '" ) = ^£glW —(N-j-ty.- < 10 > 

As a special case we mention m = 0, for which the partition function, 

vriV (-Nu 2 V 

^ri-^rE 1 ^- (ID 

3=0 J " 



is simply a truncated exponential. 

3. In this section we evaluate the Yang-Lee zeros of the partition function fllQl) . First, 
we consider the partition function as a polynomial in m. Fig. 1 shows results for fi = 0, 
/i = 0.5, and \i = 0.6. We have calculated the zeros for different values of n, e.g., n = 48, 
n = 96, and n = 192. The results for n = 192 are represented by the points in the figure. 
Of course, the exact location of the zeros is extremely sensitive to numerical round-off 
errors. Thus, the present results were obtained with the help of a multi-precision package 
pip . Typically, we performed our computations with 100-500 digits accuracy. 
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Figure 1: The zeros of the partition function in the complex m plane for /i = (upper), 
/i = 0.50 (middle) and fi = 0.60 (lower) for n = 192. The zeros of the discriminant of the 
cubic equation (^) are denoted by stars. 

The zeros fall on a curve and are regularly spacedQ. From the saddle point analysis of 
the partition function it is clear that the condition that the free energy of two different 
1 If the numerical accuracy is not sufficient, one typically observes that the line of zeros ends in a circle. 
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saddle point solutions coincides imposes a single condition on the complex variable m. 
This explains that, in the thermodynamic limit, the zeros are located on a curve in the 
complex plane (A similar argument has been given for the Ising model |4j].). If we increase 
the order of the polynomial by a factor of 2, we find that half of the zeros are close to 
those of the lower-order polynomial. The other zeros half are roughly mid-way between 
adjacent zeros of the lower-order polynomial. This leads us to the conclusion that the 
zeros become dense and lead to a cut in the complex m plane in the thermodynamic limit. 
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Figure 2: The zeros of the partition function in the complex fi plane for m = (upper), 
m = 0.30 (middle), and m = 1.0 (lower) for n = 192. The solid line in the upper figure 
represents the solution of the transcendental equation (|8]). The zeros of the discriminant 
of the cubic equation (0) are denoted by stars. Note that the scale on the x-axis of the 
lower figure is different. 

The stars in Fig. 1 represent the points at which the discriminant (|5|) of the cubic 

saddle point equation vanishes. These points coincide with the endpoints of the line of 

zeros. We have verified that the line of zeros coincides with the line at which the partition 
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function is dominated by a different solution of the third-order equation. A schematic 



picture of this line is shown in |p2]j . At the critical value of /i, the line of zeros splits into 
two lines. 

In Fig. 2 we show the zeros of the partition function in the complex /i plane for n=192 
and masses m = 0, m = 0.3 and m = 1.0. We have evaluated the zeros for n = 48 and 
n = 96 as well. Also in this case we have performed our calculations with 100-500 digits 
accuracy. The zeros are regularly spaced and their density increases homogeneously with 
n. In the thermodynamic limit, we therefore expect that they join into a cut. In all cases 
the endpoint of the line of zeros ends in a zero of the discriminant (||) (denoted by a star). 
For m = the thermodynamic limit of the line of zeros is given by the solution of the 
transcendental equation (|8|) (full line in upper figure). Chiral symmetry is broken in the 
region enclosed by this curve and is restored outside. Up to corrections of order 1/n the 
zeros coincide with this curve. For m = the density of zeros approaches zero near fi = i. 
This is not surprising since at this point the phase transition changes from first order into 
second order. 

4. In conclusion, we have evaluated the zeros of the partition function in a schematic 
random matrix model of the chiral transition at non-zero chemical potential. Since the 
coefficients of these polynomials were obtained analytically, it was possible to locate the 
zeros numerically with considerable accuracy provided that all arithmetic operations were 
executed with a precision of 100 to 500 digits. (We stress that calculations with ordinary 
double or quadruple precision arithmetic are not adequate.) We have found that these 
zeros lie on one-dimensional subsets in the complex m or fi planes. All of our results 
are consistent with a mean field analysis of the random matrix partition function. In 
particular, we have obtained an analytic expression for the curve of zeros in the complex 
fx plane. 
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Figure Captions 

Fig.l. The zeros of the partition function in the complex m plane for /i = (upper), 
\x = 0.50 (middle) and /i = 0.60 (lower) for n = 192. The zeros of the discriminant of the 
cubic equation ([?[) are denoted by stars. 

Fig. 2 The zeros of the partition function in the complex fi plane for m = (upper), 
m = 0.30 (middle), and m — 1.0 (lower) for n = 192. The solid line in the upper figure 
represents the solution of the transcendental equation (|8]). The zeros of the discriminant 
of the cubic equation (0) are denoted by stars. Note that the scale on the x-axis of the 
lower figure is different. 
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